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Abstract 

Background: In 2003, Plasmodium vivax malaria has re-emerged in central eastern China including Yongcheng 
prefecture, Henan Province, where no case has been reported for eleven years. Our goals were to detect the space- 
time distribution pattern of malaria and to determine significant environmental variables contributing to malaria 
incidence in Yongcheng from 2006 to 2010, thus providing scientific basis for further optimizing current malaria 
surveillance and control programs. 

Methods: This study examined the spatial and temporal heterogeneities in the risk of malaria and the influencing 
factors on malaria incidence using geographical information system (GIS) and time series analysis. Univariate 
analysis was conducted to estimate the crude correlations between malaria incidence and environmental variables, 
such as mosquito abundance and climatic factors. Multivariate analysis was implemented to construct predictive 
models to explore the principal environmental determinants on malaria epidemic using a Generalized Estimating 
Equation (GEE) approach. 

Results: Annual malaria incidence at town-level decreased from the north to south, and monthly incidence at 
prefecture-level demonstrated a strong seasonal pattern with a peak from July to November. Yearly malaria 
incidence had a visual spatial association with yearly average temperature. Moreover, the best-fit temporal model 
(model 2) (QIC = 16.934, P<0.001, R 2 = 0.818) indicated that significant factors contributing to malaria incidence 
were maximum temperature at one month lag, average humidity at one month lag, and malaria incidence of the 
previous month. 

Conclusions: Findings supported the effects of environment factors on malaria incidence and indicated that 
malaria control targets should vary with intensity of malaria incidence, with more public resource allocated to 
control the source of infections instead of large scale An. sinensis control when malaria incidence was at a low level, 
which would benefit for optimizing the malaria surveillance project in China and some other countries with 
unstable or low malaria transmission. 
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Background 

Malaria is one of the major causes of morbidity and 
mortality in the world, with less than one million deaths 
annually reported by WHO [1]. Malaria is a serious glo- 
bal public health problem, and its prevention and control 
is addressed in the United Nations (UN) Millennium 
Development Goals (MDG) [2]. Yongcheng prefecture, 
Henan province is a research center of the 2010-2020 
plan for national malaria elimination in China [3], which 
is an important part of global malaria action plan 
(GMAP) [4]. In the past, malaria in Yongcheng pre- 
fecture were severe, with incidence as high as 3.34 per 
100,000 population in 1970. Malaria incidence was dras- 
tically reduced after national comprehensive interven- 
tions, which included case management and Anopheles 
elimination. In 2003, malaria re-emerged in Yongcheng 
prefecture after a period of eleven years without a reported 
case [5,6]. Malaria incidence increased to 0.02 per 100,000 
population in 2006, which accounted for 4.52% of the 
total malaria cases in China [7] . 

In central eastern China including Yongcheng prefec- 
ture, where malaria incidence remains seasonal and un- 
stable, there is a need for timely confirmation of 
significant factors and development of factors targeting 
malaria interventions to curtail malaria incidence. Al- 
though studies have described malaria and its interven- 
tions in detail or factor analysis in central China [5,6,8- 
10], but thus far, attempts to develop predictive models 
of malaria epidemics with environmental variables, 
which are accurate on the local scale, have not met with 
success. 

Malaria is one of the important environmental dis- 
eases. When the environmental parameters (such as 
temperature, humidity) permit, Anopheles mosquitoes 
would transmit the pathogen, Plasmodium spp [11]. 
Therefore, to figure out how malaria varies due to sea- 
sonal or year-to-year changes in environmental variables 
is essential for the national malaria elimination plan in 
order for it to allow interventions to be adapted to the 
specific sites or times of year. Although previous studies 
proved that weather variables have exerted great influ- 
ences on malaria incidence in different regions of the 
world [12-20], controversial issues still remain. First of 
all, some experts argued that these correlations were 
questionable [20,21]. For example, some of them claimed 
that rainfall had an significant effect on the incidence of 
malaria [22-23], whereas others did not detect a signifi- 
cant relationship [24,25]. Moreover, the meteorological 
factors that have significant statistical correlations with 
malaria vary greatly between geographic areas within the 
world, which would complicates the decision-making 
process in choosing weather monitoring targets [26]. In 
addition to those controversies, Anopheles density is 
considered as a proximate factor of malaria transmission 



while climatic indicators are distal risk factors, how- 
ever, few researchers have treated mosquitoes abun- 
dance as an independent variable in statistical models 
to account for malaria incidence due to limited avail- 
ability of entomological data [27]. Last but not the 
least, studies only conducted correlation analyses might 
have been insufficient [24,28] In sum, it is necessary to 
perform further studies to elucidate the impacts of en- 
vironmental factors, such as meteorological factors, 
Anopheles density on malaria incidence; this study did 
this based on data from Yongcheng prefecture, China 
as one example. 

The data used in this study were collected from Yong- 
cheng prefecture during the period between 2006 and 
2010, and were analyzed using Geographic Information 
System (GIS) and Generalized Estimating Equations 
(GEE) approach. The aim was to provide not only a sci- 
entific basis for malaria monitoring and control in Yong- 
cheng prefecture but also valuable information for 
malaria elimination in other areas of seasonal and un- 
stable malaria incidence. 

Methods 

Study area 

The site selected for the study was Yongcheng prefec- 
ture, located in eastern Henan Province (latitude 33° 
42'~34°18' N and longitude 115°58' ~ 116°39' E) 
(Figure 1). The prefecture has an area of 1994.49 km 2 
and a population of 1.46 million. Yongcheng prefecture 
has a temperate and monsoonal climate with four clearly 
distinct seasons, with monthly maximum air temperate 
of 27.3°C, monthly minimum temperature of -1.9°C, 
and monthly average precipitation of 219.05 mm. The 
geographic landscape and complex climate situation 
such as suitable temperature and humidity, abundant 
rainfall, and existence of water bodies, provided 
favourable breeding sites for Anopheles, which would 
contribute to making it a suitable environment for P. 
vivax malaria transmission. 

Data collection and management 

In China, malaria is a statutory notifiable category B in- 
fectious disease [30], reported to the National Chinese 
Center for Disease Control and Prevention (CDC). 
Monthly malaria case data from 2006 to 2010 were 
obtained from Yongcheng CDC. Physicians at hospitals 
and clinics reported suspected malaria cases to the local 
CDC, and the cases were defined based on the diagnos- 
tic criteria and principles of management for malaria 
(GB 15989-1995) issued by Ministry of Health of the 
Peoples Republic of China. Only the cases confirmed 
clinically and by laboratory test, including thick and thin 
blood smear, were included in our study. The population 
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data were obtained from reports produced by Yong- 
cheng Statistical Bureau. 

Monthly weather data from 2006 to 2010 were 
obtained from China Meteorological Data Sharing Ser- 
vice System (http://cdc.cma.gov.cn/), including monthly 
average temperature (T avg ), maximum temperature (T max ), 
minimum temperature(T min ), average relative humidity 



(H avg ), minimum humidity(H min ), total 24-hour rainfall 
(R), average wind velocity(W avg ), maximum wind velocity 
(W max ), extreme wind velocity(W e ) and duration of 
sunshine (S). Due to the unavailability of weather data 
of Yongcheng prefecture, geostatistical methods of 
kriging was used to estimate values of meteorological 
parameters in Yongcheng using monthly weather 



Zhang et al. BMC Public Health 2012, 12:544 
http://www.biomedcentral.com/1471-2458/12/544 



Page 4 of 13 



observations measured at its eight neighbouring climatic 
stations, including Shangqiu, Xihua, Dangshan, Xuzhou, 
Haozhou, Suzhou, Fuyang, and Bengbu weather station. 
The distribution of those weather stations was mapped 
in Figure 2. 

Data on Anopheles density were collected by Yong- 
cheng CDC. Since 2005, when Maqiao village, Yong- 
cheng prefecture was identified as a national malaria 
monitoring site, samples of Anopheles have been col- 
lected fortnightly from June to August using China mal- 
aria surveillance protocol [30]. Anopheles density was 
measured by two methods: human-bait catches and bed- 
net trap. Human baits were positioned in outdoor areas 
near households and mosquito breeding sites to collect 
mosquitoes continuously from 7:00 PM to 10:00 PM. 50 
bed-net traps were distributed to 50 households to col- 
lect indoor mosquitoes from 6:00 AM to 7:00 AM. 
Afterwards, the collected mosquitoes were then sepa- 
rated, sorted, and identified. All collected Anopheles in 
the study period were morphologically identified as An. 



Sinensis, and a sample was re-identified as An. sinensis 
by PCR, which is in accordance with the results of previ- 
ous studies that An. sinensis was the only species 
reported to be the vector capable of transmitting P. 
vivax in Yongcheng [5,8]. Semi-monthly An. sinensis 
density measured by human-bait catches (D bait ) and 
bed-net trap (D net ) was calculated as the number of 
mosquitoes per person*night or per 50 mosquito nets. 
The monthly mean value for An. sinensis density used in 
this study was calculated from semi-monthly data. If 
monthly An. sinensis density is unavailable in some 
month (except for June to August) of the year, we trea- 
ted it as zero in the process of figure drawing and model 
building because of its seasonality. 

Ethical approval 

Ethical approval of this study was obtained from the 
Ethical Committee of China CDC and permission was 
also got from the Municipal Government, the Municipal 
Health Bureau and CDC in Yongcheng city. 



Xihuan 



Legend 

O Meteorological site 
| Yongcheng prefecture 
the neighbourhood of Yongcheng 

NAME 

] Anhui province 
] Henan province 
] Jiangsu province 
| Shandong province 




Figure 2 Positions of the eight climatic stations neighbouring Yongcheng prefecture. 



Zhang et al. BMC Public Health 2012, 12:544 
http://www.biomedcentral.com/1471-2458/12/544 



Page 5 of 13 



Kriging analysis for spatial interpolation of 
meteorological factors 

Some weather values at unsampled locations were inter- 
polated by Kriging method [31]. Monthly observed 
values of climatic data from eight sampled weather sta- 
tions neighbouring Yongcheng prefecture was the main 
input variables. Kriging was performed following the 
procedures described by Journel et al. [32] and Burgess 
et al. [33], which included preliminary data analysis, struc- 
tural data analysis, log kriging estimations, and image 
generations of spatial results. Spatially distributed values 
of weather factors in Yongcheng prefecture were esti- 
mated based on spherical model in this study, which 
was regarded as the most widely used semi-variogram 
model [34]. Monthly weather values at county- or 
prefecture-lever were calculated by using the method of 
Zonal statistics. The above analyses were conducted with 
ArcGIS 9.2 software (ESRI Inc, Redlands, California). 

Mapping malaria incidence and average temperature with 
GIS 

We conducted GIS -based analyses of the spatial dis- 
tribution of yearly malaria incidence and average 
temperature as well as their visual correlation. The 
annualized average incidence of malaria per 100,000 
persons and annual average temperature at each tow- 
lever over the five years (from 2006 to 2010) were cal- 
culated and mapped based a town-level polygon map 
of Yongcheng at 1:1000 000 scale with ArcGIS 9.2 
software (ESRI Inc, Redlands, California). Regions with 
different intensities of malaria incidence and various 
values of average temperature were marked with differ- 
ent colours on the town-level, with higher incidence 
and temperature being indicated by darker colour. 

Temporal analysis with GEE technology 

The monthly malaria incidence, An. sinensi density, and 
weather variables were calculated and plotted to observe 
their seasonal fluctuations and correlations from 2006 to 
2010. 

Time varying influencing factors were treated with dif- 
ferent time lags, from 0- to 3- month lags, to account 
for delays in their effects on malaria incidence. The lag 
size was determined by comparing quasi-likelihood 
under the independence model criterion (QIC) values in 
models with various lag sizes [35]. 

Univariate analyses were made by regressing single 
factors of interest against monthly malaria incidence to 
estimate crude associations between malaria incidences 
and influencing factors. Multivariable models were built 
to examine the effects of combinations of influencing 
factors on malaria incidence. Candidate factors selected 
for the multivariable model were determined through 



statistical performance of factors in the univariate ana- 
lysis, and hypothesized relationships. Candidate influen- 
cing factors were inputted into the model in their 
presumed order of importance, and then non- 
statistically significant factors were removed in their pre- 
sumed inverse order of importance unless the remaining 
factors were deemed important for theoretical reasons at 
a = 0.05 level. 

Two kinds of multivariable models were built up in 
this study using GEE approach. The difference between 
two models lay in whether they included malaria inci- 
dence of the previous month into the model. Model 1 
was constructed to estimate the relationships between 
An. sinensis density, weather variables and malaria inci- 
dence. Model 2 was developed to examine the effects of 
An. sinensis density, weather variables and the malaria 
incidence of the previous month on malaria incidence. 

The goodness of fit of the GEE model was measured 
by "marginal R-square", which was interpreted as the 
amount of variance in the response variables that were 
explained by the fitted model [36], and QIC value, which 
was useful in selecting an appropriate correlation struc- 
ture [8,37,38]. The model with a lowest QIC score and a 
highest R-square was preferred. GEE analysis was imple- 
mented by STATA software 11.0 (Stata Corp. College 
Station, Texas). 

Results 

Spatio-temporal distribution of malaria incidence and 
influencing factors 

From 2006 to 2010, 6,546 malaria cases were reported 
(7.88 cases per 100,000 population) in Yongcheng prefec- 
ture, which were all diagnosed as P. vivax malaria. The 
majority of cases were farmers (76.08%). The annual mal- 
aria incidence at the town level showed an increasing 
pattern from the north to south of Yongcheng prefecture 
between 2006 and 2010 (Figure 3). The areas of highest 
malaria incidence were distributed in the southern region 
of Yongcheng prefecture, Maqiao village and Lizhai vil- 
lage, bordering northern Anhui Province. As Figure 4 
showed, the yearly average temperature also rose from 
the north to south at the town level, with higher values 
distributed in the Southeast of Yongcheng prefecture. 

There was a decline (15.45 cases per 100,000 popula- 
tion) in annual malaria incidence from 2006 to 2010, 
and the highest incidence occurred in 2006 with 16.72 
cases per 100,000 population. Monthly P. vivax malaria 
incidence showed a seasonal pattern, whose peak period 
was from July to November, a period when nearly 
86.58% of total malaria cases were reported (Figure 5). 
The vector specie that transmitted malaria in Yongcheng 
prefecture was morphologically identified as An. sinensis, 
with PCR results of a random sample of 30 mosquitoes 
confirmed as An. sinensis. An. sinensis density measured 
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by different methods particularly the human-bait collec- 
tions displayed a seasonal variation where the valleys 
and peaks closely corresponded to malaria incidence. 
Moreover, meteorological factors also fluctuated season- 
ally during the study period in addition to precipitation, 
showing possible relationships with malaria incidence 
(Figure 5). 

Crude correlations between malaria incidence and 
environmental factors 

Univariate analysis was preformed to provide a basis for 
selecting significant factors for malaria incidence from 
twelve influencing factors. Monthly An. sinensis density, 



temperature, humidity, and rainfall had significant posi- 
tive correlations with malaria incidence with delays of 
zero to three months, while wind velocity showed nega- 
tive correlations with delays of zero and one month. 
There was not a significant association observed be- 
tween duration of sunshine and malaria incidence 
(Table 1). All these significant factors with various lag 
sizes were entered into the multivariable model. 

Model building and evaluation 

Twelve variables, such as D bait , D net , T avg , T max , T min , R, 
H avg > H min , W avg , W max> W e , and S were inputted into a 
model for variable selection. The model with n-binomial 
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Figure 4 Yearly average temperature at county level from 2006 to 2010 in Yongcheng prefecture, China. The county-lever areas are 
colour-coded according to yearly average temperature, with higher temperature distributed in the Southeast of Yongcheng prefecture. 



distribution (model 1) effectively explained the relation- 
ships between An. sinensis density, weather variables and 
malaria incidence (Table 2, model fitting QIC = 37.141, 
R 2 = 0.520). Significant factors contributing to malaria 
incidence showed in model 1 were D bait (p = 0.050, 
P = 0.033) at 1 -month lag, T max at 1 -month lag 
(p = 0.174, P<0.001), and W e at 0-month lag (p = -0.145, 
P = 0.009) (Table 3). No interaction terms were signifi- 
cant in the multivariable model. 
Meanwhile, thirteen variables, including D bait , D net , 



T T 



:> T m j n , H av g, H n 



R, W av2 , W max , W e , S and 



the malaria incidence of the previous month (I p ) were 
added into another model to select variables. A best-fit 
model (model 2) was with n-binomial distribution 
(Table 2, model fitting QIC = 16.934, R 2 = 0.818), show- 
ing that T max at a lag of one month (p = 0.153, P<0.001), 
H avg at one month lag (p = 0.024, P = 0.010), and I p 



(p = 0.049, P< 0.001) significantly associated with malaria 
incidence (Table 4). However, D bait and D net failed to be 
included in model 2. No interaction terms were signifi- 
cant in the multivariable model. 

Expected values fitted by model 1 or model 2 and ac- 
tual observed values of malaria incidence in 2006-2010 
were shown in Figures 6 and 7. Fitted values predicted 
by model 1 were in agreement with actual values of mal- 
aria incidence in the first three years (F = 0.645, 
P = 0.425), but they were over-estimated to a greater ex- 
tent in the last two years (F = 29.052, P<0.001). How- 
ever, fitted values predicted by model 2 were in 
accordance with actual values of malaria incidence in 
the study period of 2006-2010(F = 0.032, P = 0.857). 

Scatterplots of the differences between the predicted 
and observed malaria incidence were plotted (Figures 8 
and 9), and the arithmetic mean and 95% confidence 
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Figure 5 Time series plots of malaria incidence and influencing factors from 2006 to 2010 in Yongcheng prefecture, Henan. Influencing 

factors include D bait; D net; T avg , T max; 

Tmin/ H av g, H m j n , 

R, w avg , w max , w e . 



interval for the differences were calculated. The average 
difference values between actual values and fitted values 
of model 1 (D-value 1) or model2 (D-value 2) were re- 
spectively 1.476/100,000 (9596CI: -1.498, 4.451/100,000), 
0.769/100,000 (9596CI: -1.347, 2.885/100,000). There was 
no significant difference between the actual values and 
the values predicted by model 1 or model 2. 

The best-fit model was model 2, because it had not 
only a smaller D-value (0.769/100,000), but also a smal- 
ler QIC value and a larger R 2 (QIC = 16.934, R 2 = 0.818), 
compared to model 1 (D-value = 1.476/100,000, QIC = 
37.141, R 2 = 0.520). 

Discussion 

This is the first study spatially and temporally exploring 
the effect of weather variables and vector parameters on 
malaria incidence in China. The data involved in this 
study were relevant because all of them were obtained 



from national monitoring data. In this study, clear 
spatial heterogeneity and temporal clustering of malaria 
incidence could be found, with higher incidence distrib- 
uted in the Southern Yongcheng spatially and in July to 
November temporally. The finding indicated that areas 
and months with higher malaria transmission risk 
should be focused on more public health attention and 
resources. Spatial heterogeneity of malaria incidence 
(higher incidence in the south) could be explained by 
the spatial variability of temperature and the distribution 
of malaria imported cases. First of all, Southern Yong- 
cheng prefecture has been confronted with higher pres- 
sure from imported cases, as its was adjacent to Suixi, 
Guoyang, Xiao, and Huaibei county in Anhui province, 
which were identified as high-endemic areas of malaria 
(incidence > 30/100,000) [10]. Furthermore, its spatial 
decreasing pattern from the north to south was some- 
how in accordance of temperature change, particularly 
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Table 1 Coefficients of univariate analysis among environmental factors on P. vivax malaria incidence 



Lag-time 
(months) 


Distribution 
type 


Dbait 


Dnet 


"l"avg 


Tmax 


Tmin 


Havg 


Hmin 


R 


Wavg 


^max 


W e 


S 


0 


Poisson 


0.091*** 


0.164 


0.074*** 


0.078*** 


0.073*** 


0.073**** 


0.044*** 


0.002 


-2.782*** 


-0.453** 


-0.201 


-0.005 




n-binomial 


0.144*** 


0.186 


0.122*** 


0.113*** 


0.110*** 


0.107** 


0.068*** 


0.002 


-3.180*** 


-0.391* 


-0.167** 


-0.010 


1 


Poisson 


0.120*** 


0.250* 


0.143*** 


0.148*** 


0.137*** 


0.097*** 


0.059*** 


0.006** 


-1.782 


-0.185 


-0.040 


-0.003 




n-binomial 


0.222*** 


0.370** 


0.192*** 


0.199*** 


0.188*** 


0.110*** 


0.073*** 


0.008** 


-2.073*** 


-0.150 


-0.035 


-0.004 


2 


Poisson 


0.116*** 


0.248* 


0.184*** 


0.188*** 


0.167*** 


0.086*** 


0.064*** 


0.009*** 


-0.320 


0.121 


0.109 


0.003 




n-binomial 


0.173*** 


0.398*** 


0.181*** 


0.185*** 


0.175*** 


0.067** 


0.049 


0.015*** 


-0.404 


0.126 


0.128 


0.002 


3 


Poisson 


0.081** 


0.079 


0.124*** 


0.130*** 


0.111*** 


0.039 


0.039 


0.007*** 


0.950* 


0.469*** 


0.199*** 


0.009 




n-binomial 


0.105** 


0.130 


0.128*** 


0.130*** 


0.124*** 


0.032 


0.023 


0.012*** 


1 .392*** 


0.324** 


0.277*** 


0.007 



("*": The regression coefficient is significant at level of a = 0.05, "**": The regression coefficient is significant at level of a = 0.01, "*' : 
significant at level of a = 0.001). 



*": The regression coefficient is 



in 2006. Only temperature was selected to assess the in- 
fluencing factor associated with spatial heterogeneity 
owing to the founding from previous studies that 
temperature might be a major determinant of malaria 
incidence in China [26]. Therefore, public resource 
should be allocated proportionally in different areas 
based on the risk predicted by imported possibility and 
average temperature, and close collaboration should be 
established between Henan and Anhui to effectively con- 
trol and prevent malaria together. From the perspective 
of time, annual malaria incidence in Yongcheng prefec- 
ture showed an obvious decrease from 2006 to 2010. 
However, prevention and control of malaria should not 
be taken slightly, because Yongcheng prefecture is still at 
risk for malaria due to the existing of An. sinensis, suit- 
able weather condition for the growing of An. sinensis 
and Plasmodium, and the risk has been possibly varying 
with weather changes and population movement [39,40]. 
Determining the principal influencing factors of malaria 
incidence would be beneficial for malaria risk assessment 
and thus providing a basis for the policy making for mal- 
aria control technologies. 

Malaria is a vector-borne infectious disease, and as such, 
is sensitive to environmental change [40-42]. Anopheles 
density, as a proximate environmental factor of malaria 
transmission, play an important role in estimating and 
predicting malaria risk [27]. Climatic variables have also 
been established as important environmental drivers of 
malaria transmission [43], because of their impacts on 
the growth and reproduction rates of mosquitoes, the 
temporal activity pattern of the population as well as the 
life cycle of Plasmodium [44-47]. 

Table 2 Comparison of GEE models 





Model 1 


Model 2 




QIC 


QIC_u 


QIC 


QIC_u 


Poisson Distribution 


131.941 


119.726 


351.132 


319.353 


n-binomial Distribution 


16.934 


23.029 


37.141 


41.236 



The best-fit model (model 2) derived from the study 
was reliable and had a good fit and predictive validity 
(QIC =16.934, P<0.001, R 2 = 0.818), which provided 
insights into the most important drivers of P. vivax mal- 
aria, including maximum temperature, average humidity 
and incidence of previous month that influenced sea- 
sonal fluctuation of P. vivax malaria incidence. The re- 
sult that temperature rise would contribute to malaria 
transmission was in agreement with some researchers 
[46-52], although there were still some other researchers 
who argued that this relationship was not significant 
[19,53], or that it was uncertain [54]. It has been demon- 
strated that temperature increase would improve the 
survival chances of Anopheles and thus contribute to the 
malaria transmission [55,56] Moreover, relative humidity 
exerted an influence on the survival of mosquito eggs 
and adults and the moderate increase in malaria risk 
associated with average humidity observed in this study 
was consistent with previous findings [57]. Conversely, 
some literature found a correlation between rainfall and 
malaria [49], while other studies found no correlation 
[51]. In this study, rainfall failed to enter the best-fit 
model as a predictor of malaria epidemics. The 
phenomenon could be explained by the complex non- 
linear association between rainfall and malaria incidence. 
Rainfall is beneficial to the growth and reproduction of 

Table 3 Poisson and n-binomial regression of malaria 
incidence on influencing factors (the result of Model 1) 

Lag-time Variables Poisson Distribution n-binomial Distribution 



(months) 


P 


SE 


P 


P 


SE 


P 


1 D bait 


0.038 


0.017 


0.025 


0.05 


0.023 


0.033 


1 T max 


0.096 


0.017 


<0.001 


0.174 


0.015 


<0.001 


1 H avg 


0.044 


0.011 


<0.001 








0 W e 








-0.145 


0.055 


0.009 


constant 


-3.811 


0.958 


<0.001 


-0.975 


0.726 


0.179 



{Coef: regression coefficient, SE: estimated regression coefficient standard 
error). 
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Table 4 Poisson and n-binomial regression of malaria 
incidence on influencing factors (the result of Model 2) 

Lag-time Variables Poisson Distribution n-binomial Distribution 



(months) 




P 


SE 


P 


p 


SE 


P 


1 


Tmax 


0.137 


0.019 


<0.001 


0.153 


0.012 


<0.001 


1 


Havg 








0.024 


0.009 


0.01 


2 


R 


0.005 


0.001 


<0.001 








0 




0.039 


0.003 


<0.001 


0.049 


0.006 


<0.001 




constant 


-2.51 


0.488 


<0.001 


-4.449 


0.589 


<0.001 



{Coef: regression coefficient, SE: estimated regression coefficient standard 
error). 



mosquito if it is moderate, because it often leads to pud- 
dles and increased local humidity; however, excessive 
rain can also wash away eggs and completely destroy 
breeding sites [58]. This result indicated that it was not 
necessary to consider rainfall as a predictor in Yong- 
cheng, which made malaria surveillance simpler in this 
area. 

In addition to factors mentioned above, although An. 
sinensis densities failed to enter into best-fit model 
(model 2), D bait was included in model 1. The reasons 
why D bait rather than D net were included in model 1 
probably lay in two aspects. Firstly, An. sinensis is 
slightly exophagic (biting outdoors) [59], and thus D bait 
is more representative of malaria transmission in Yong- 
cheng. Secondly, most of the malaria cases in Yongcheng 
were farmers, and they would like to sleep and work 
outdoors in the summer without effective protections, 
thus having more opportunities to be infected [5,9]. 



Therefore, D bait would contribute more to the prediction 
of malaria incidence. Moreover, comparing the time 
series plots of D bait (Figure 5), expected values predicted 
by model 1 (Figure 6) and actual values (Figure 6), we 
found out expected values were over-estimated in 2009 
and 2010 when actual malaria incidence was relatively 
low, which mainly due to the existence of predictor D bait 
reasoned by their similar rising tendency. It could be 
concluded that An. sinensis density probably had its 
shortcomings as a routine monitoring and predicting 
index. An. sinensis density may be a better predictor of 
malaria incidence when transmission is relatively high, 
as many of the female mosquitoes may have been 
infected by Plasmodium, and increase in An. sinensis 
density would lead to a direct rise in malaria incidence. 
However, when the incidence is at a low level, most of 
the female mosquitoes are possibly free of Plasmodium. 
In this circumstance, when a healthy human being is bit- 
ten by a female An. sinensis the probability of infection 
by Plasmodium is low. Therefore, we concluded that rise 
in An. sinensis density would possibly contribute little to 
the increase of malaria incidence in low transmission 
areas, which agreed with the result of another study that 
malaria transmission potential would be very low in 
spite of a high human biting rate in unstable malaria 
areas [60] . Therefore, control target should vary with the 
severity of malaria epidemics. When malaria incidence is 
high, public resource allocation should be focused on 
mosquito control and elimination; however, when mal- 
aria incidence is low, the key control point should lie in 
controlling sources of infections. Furthermore, it is 
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Figure 6 Expected values of model 1 and observations for malaria incidence in the study period of 2006-2010. 
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Figure 7 Expected values of model 2 and observations for malaria incidence in the study period of 2006-2010. 



necessary to find a substitute, such as entomological in- 
fection rate (EIR), which can overcome the weakness of 
An. sinensis density as an indicator for malaria surveil- 
lance and prediction, although it can be used as a good 
index for predicting malaria potential risk as a previous 
study showed [40]. 

As far as the lag effect was concerned, this study found 
significant one month lag effects of entomological and 
meteorological variables on malaria incidence, and this 
finding was supported by several earlier studies 
[17,19,35,49]. This phenomenon could be explained by 
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Figure 8 Scatterplot of difference between expected values 


predicted by model 1 and observed values of malaria incidence 


in the study period of 2006-2010. 



the approximately one month duration of malaria infec- 
tion cycle. The time incorporates several processes 
above. An adult mosquito first bites an infected human, 
and then the parasite develops in the adult mosquito 
(Extrinsic Incubation Period). Ten days later when the 
P. vivax sporozoites move the salivary glands, the mos- 
quito transmits malaria to a human when it takes another 
blood meal. Once the person is infected, time to devel- 
opment of malaria symptoms and infectivity (Intrinsic 
Incubation Period) takes about another 1-2 weeks 
[40,61]. Knowing the approximate lag size of effects on 
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Figure 9 Scatterplot of difference between expected values 


predicted by model 2 and observed values of malaria incidence 


in the study period of 2006-2010. 
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malaria incidence would benefit us to get prepared for 
quick and effective response on malaria epidemic easily 
at least one month in advance. 

The transmission of malaria is complicated, and we 
still need further research to figure it out. For example, 
the temporal variation in malaria incidence could be also 
partially explained by continuous and effective control 
efforts Chinese local and national health agencies made, 
such as treatment in the rest period of malaria(con- 
ducted from 2004) [62] and comprehensive vector con- 
trol action characterized by biological larviciding and 
residual spray(conducted from 2007) [8]. Human inter- 
ventions are failed to be inputted into the predictive 
model in this study because it is difficult to measured 
and quantify. 

Conclusions 

Areas and months with higher malaria transmission risk 
should be focused on more public health attention and 
resources. The model developed in this study success- 
fully predicted the expected incidences of malaria based 
on historical malaria epidemics and a combination of 
weather factors at one month lag, which would simplify 
malaria surveillance by targeting control of malaria more 
effectively. Furthermore, we concluded that (l)more ef- 
fective indicators, such as weather variables and Plasmo- 
dium infection rate of mosquitoes should be considered 
to further optimize current malaria monitoring and con- 
trol methods; (2) and malaria control targets should vary 
with intensity of malaria incidence, with more public re- 
source allocated to control the source of infections in- 
stead of large scale An. sinensis control when malaria 
incidence was at a low level, which would benefit for op- 
timizing the malaria surveillance project in China, and 
also could be used for the malaria monitoring and early- 
warning in some other countries with unstable or low 
malaria transmission. 
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